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We explore a long Josephson contact transporting Cooper pairs between ID charge-neutral chiral 
Majorana modes in the leads via charged Dirac chiral modes in the normal region. We investigate 
the regimes of (i) transparent contacts and (ii) tunnel junctions implemented in 3D topological 
insulator/superconductor/magnet hybrid structures. The setup acts as a SQUID controlled by the 
magnetic flux enclosed by the chiral loop of the normal region. This chirality leads to the fractional 
h/e-periodic pattern of critical current. The current-phase relation can have sawtooth-like shape 
with spikes at unusual even phases of 27rn. 


I. INTRODUCTION 

Intensive studies of Majorana fermion physics in solid 
state materials were motivated by possible applications 
of these states in topological quantum computation^’^. In 
condensed matter, initially, Majorana fermions were pre¬ 
dicted to exist in spinless p-wave superconductors (SC) 
and fractional QHE^“^. Later, the realizations of Ma¬ 
jorana bound states (MBS)^“^ and ID charge-neutral 
modes^ were proposed in hybrid structures where heli¬ 
cal electronic system has a proximity effect with exter¬ 
nal s-wave pairing potential and Zeeman exchange field. 
The helical systems with spin and momentum locking 
are realized in surfaces and edges of topological insu¬ 
lators (TI), quantum spin-Hall films^^ and semiconduc¬ 
tors with strong spin-orbital coupling^’^^. The transport 
signatures of MBS are associated with anomalous Jtt- 
periodic Josephson effect^’^’^^’^^. On the other hand, 
charge-neutral chiral Majorana modes (yMM) which are, 
for instance, supported by SC-magnet junctions on 3D 
TI surface^’^^, are expected to show unconventional in¬ 
terferometry if involved in charge transfer between nor¬ 
mal metal leads. Various yMM-based devices operating 
as Mach-Zehnder^^’^^, Fabry-Perot^^’^^, and Hanbury- 
Brown Twiss^^ quantum interferometers have been pro¬ 
posed. In context of quantum computation they can 
serve as a readout tool of qubit states encoded by MBS 
in vortex cores^^. 

In the paper we calculate DC Josephson current-phase 
relation (C^R) of a long spinless contact between ID 
gapless yMMs. The system under consideration repre¬ 
sents a chiral Fabry-Perot interferometer implemented on 
a surface of 3D TI partially gapped by SC and magnetic 
(M) islands. Similar ideas were explored in Refs. 14- 
18. In these works normal dissipative transport of Dirac 
fermions influenced by their splitting to Majorana modes 
and the interference of the latter were studied. In con¬ 
trast, we study the equilibrium phenomenon of Joseph¬ 


son current carried by the Andreev states formed in the 
chiral Dirac liquid. The latter connects between two su¬ 
perconducting leads, which are effectively gapless due to 
chiral Majorana modes at their edges. Spinless structure 
of the junction follows from the absence of spin degen¬ 
eracy of 2D helical surface states. Presence of Zeeman 
field breaks time-reversal symmetry and produces gap¬ 
less chiral channels at the magnetic domain boundaries. 
In our system the normal conducting region consists of 
two separated chiral Dirac channels supported by mag¬ 
netic domain walls. 

We explore two models of (i) transparent and (ii) 
tunnel junctions. The first one invoives four chirai Y- 
junctions, operating as direct charged-to-neutrai fermion 
converters. The second one has two tunnei contacts 
formed by thin layers of magnetic material where Ma¬ 
jorana and normal Dirac channels are hybridized. The 
length of the N-region, formed by two counter propa¬ 
gating ID Dirac modes, is assumed to be larger than 
coherence length of the induced superconductivity. In 
other words, the Thouless energy of the normal conduct¬ 
ing part of the system under consideration is significantly 
smaller than the proximity gap. 

We take into account a contribution to the transport 
from sub-gap ID chiral states only, neglecting 2D bands. 
Since our superconducting leads are gapless due to the 
presence of yMMs, the spectral current is continuous and 
consists of smeared Andreev levels. Assuming that su¬ 
perconducting leads are large and have a fixed chemicai 
potentiai, we obtain the conventionai 27r-periodicity of 
the C4>R. In other words the non-equiiibrium 47r-periodic 
Josephson effect (which is predicted^^ for systems with 
zero-energy MBSs) turns out to be irrelevant here. 

The chiral structure of the normal part consisting of 
two ID spinless Dirac modes allows to make Andreev 
pairs non-local. This means that the pair resides on 
two separated Dirac channels and one can apply mag¬ 
netic flux / inducing Aharonov-Bohm phase (^ab = 
27Tfl{hle) = 7r//4>o, where 4>o = hjifle). It follows 
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that the critical current has a fractional 2^o-periodicity 
rather than the conventional one. The realization of 
the normal region by two separated chiral channels in 
the present setup is one of the most relevant distinctions 
from the junctions with 2D topological insulator where 
Andreev pairs exist on the same edge^^’^^. 

The paper is organized as follows. In Section II we 
present main ideas of yMMs realization in hybrid struc¬ 
tures on a surface of 3D topological insulator. In Section 
III we consider the regime of transparent Josephson con¬ 
tacts . We derive the S'-matrix of the contacts by means 
of matching the field operators in the chiral channels. 
Next we calculate Josephson C^R at arbitrary tempera¬ 
ture and flux enclosed by Dirac modes. The Section IV 
addresses the tunnelling regime, in which Major ana and 
Dirac modes meet at the tunnel contacts formed by mag¬ 
netic constrictions. In Section V we discuss the Majorana 
representation of charged-to-neutral fermion scattering in 
transparent and tunnel contacts. Section VI summarizes 
the results of this paper. The technical details are pro¬ 
vided in Appendices A and B. 


II. MAJORANA MODES IN 3D TI HYBRID 
STRUCTURE 

As was proposed by Fu and Kane^, the yMM is sup¬ 
ported by SC/magnet/3D-TI structure like that shown 
in Figure 1, with an effective low-energy ID Hamiltonian 
following from a solution of a 2D Bogolyubov-de Gennes 
equation. The 2D surface states are described by the 
Rashba Hamiltonian. The system is assumed to have 
a chemical potential /i = 0, corresponding to the Dirac 
point. The induced 5-wave SC pairing potential is given 
by in the ^ < 0 half plane, while at ^ > 0 there is 

magnetic material inducing Zeeman exchange field Maz. 
The full Hamiltonian reads 

dxdy^+h^ , 

h = ivTz’L • (cr X V) — M0{y)(Jz 
+ (r+Ae“*^®^+T_Ae*^®^)6>(-2/), (1) 

where the field operator of the 2D surface is T = 
7 /;^,7/^^, — 7 /;^]^, and cr and r are Pauli matrices in 
spin and Nambu spaces, respectively. The field T sat¬ 
isfies the charge conjugation constraint T = 
while the eigenstates obey = dyTy^'^p. The low-energy 
\ep\ < A,M eigenvalues of the Bogolyubov-de Gennes 
equation h^p = Sp^p are single degenerate and correspond 
to a ID chiral mode with the linear spectrum 

Sp = sign(M)up. (2) 

Single degeneracy of eigenvalues implies that Bogolyubov 
quasiparticle operator 


( 3 ) 


represents a Majorana mode obeying Xp = X-p- The 
corresponding wave function of the yMM reads 


pdf-- 




—sign(M)e + D 
—sign(M)e*^ ^ ^ 


e 


ipx 


( 4 ) 

The momentum p here is directed along x-axis and 
g^^\y) = !^2{lsc + Im) is the trans¬ 

verse shape of ID guiding channel. The coherence lengths 
(transversal decay lengths) are given by I sc = hv j /X and 
= hvjM. The superscript (M) in Eqs. (3) and (4) 
emphasizes the Majorana nature of the mode. 



FIG. 1. Superconductor/magnetic insulator boundary on the 
surface of 3D topological insulator. The boundary supports 
a chiral Majorana mode with the chirality depending on the 
sign of the magnetization. 

Another building block of chiral interferometers is a 
domain wall on a surface of a 3D TI where the magneti¬ 
zation sign is changed. If we consider such a —M/+M 
boundary along the x-axis, which is described by the 
term Mcr;^sign(^) in the Bogolyubov-de Gennes Hamil¬ 
tonian, we end up with the Dirac chiral mode with the 
same spectrum found for yMM (2). In the Nambu no¬ 
tation any of 5 eigenvalues are doubly degenerate and 
related to orthogonal electron and hole eigenstates 
and ^p^^ = (7yTy^^_^p*, where 


ipHx,y) = 


1 

isign(M) 

0 

0 


( 5 ) 


and = exp( —1^|//^)/\/2/^. There are two indepen¬ 

dent excitations with energy 5 in Nambu notation in this 
—M/+M case - the electron with momentum p and Bo¬ 
golyubov operator ^jjp and the —p hole with This 


Xp = 


'x,y) • ^{x,y) 
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field is complex ijjp ^ and corresponds to a charged 
mode. We introduce here ID operators '0(x),x(x) 

X{x) = x+(a;) = J i;{x) = J 

with the anticommutation rules given by 

{tp+ixi), 1p{X2)} = S{xi - X2), 

{lp{xi),'lp{x2)} = 0 , 


III. TRANSPARENT REGIME 
A. S'-matrix of a transparent contact 

In this section we consider the Josephson junction 
shown in Figure 2. This consists of two counterpropagat- 
ing chiral Dirac modes which scatter at the left and right 
contacts with the superconducting leads. Each of these 
contacts consists of a pair of chiral Y-junctions where 
electrons convert into a pair of Majorana fermions with 
opposite chiralities or vice versa (see Figure 3). The full 


{x(a:i), X{x2)} = S{xi - X 2 ) ■ 

The same relations hold in the Heisenberg picture at 
equal times ti = ^ 2 - After that the secondary quantized 
operators in the low energy range jej < A, M can be 
found as sums over subgap chiral states. For the charge- 
neutral Majorana excitations we obtain 



^M(x,y) = j ^^l^\x,y)xp, 

while for the Dirac mode the field reads 

'^D{x,y) = y ^ {^pHx,y)ijp + ^^p^\x,y)i;1:p^ . 

Substituting these fields into the second quantized 
Bogolyubov-de Gennes Hamiltonian and performing the 
transverse integration we obtain the following ID Hamil¬ 
tonians: 

iv f 

Hm = sign(M)— / dxx{x)d^x{x) (6) 

and 

H]j = sign{M)iv J dx^^{x)dxi’{x). (7) 

The Hamiltonians (6) and (7) correspond to coherent 
propagation of the excitations in ID guiding channels 
with the Fermi velocity v and chirality dependent on 
the sign of the magnetization. The 1/2 in the Majo¬ 
rana Hamiltonian Hm reflects the fact that independent 
excitations in the yMM can be considered either at pos¬ 
itive or negative energies only. Say, the bottom branch 
of the chiral mode (2) at p < 0 is redundant. 

The wave functions or show that the spin 

direction is locked to the momentum. More specifically, 
the spin textures of the guiding channels of Majorana 
and Dirac modes are orthogonal to the momentum direc¬ 
tion in the particular case of Rashba type Hamiltonian 
(1). As a consequence the spin textures are antiparal¬ 
lel in the counterpropagating channels. More generally, 
any scattering at a junction will be accompanied by the 
corresponding spin rotation. 


FIG. 2. Scheme of transparent Josephson junction on the 
surface of 3D topological insulator. Red lines stand for gap¬ 
less Majorana fermion channels and arrows reflect chiralities. 
Superconducting electrodes, marked as green, have phase dif¬ 
ference Light and dark grey bars are magnetic materi¬ 
als which induce exchange fields of the opposite polarizations 
±M. The line of the sign change of M supports chiral charged 
modes marked as blue. Magnetic flux / in — M region induces 
Aharonov-Bohm phase (j)AB — 7r//$o- 



'^out “hM^ ^ 

Vi 

SC ' 

^sc 

Xout \ 

1 

^ -M 

d 

1 

J bm ^ 


FIG. 3. Structure of the contact. Black arrows shows spin 
texture of the chiral modes. 

S'-matrix of this contact is derived from two operator re¬ 
lations for both upper and lower Y-j unctions, described 
by Sout and Sin matrices found in Refs. 14 and 15. Let 
us assume first that 4>5C' = 0 in the electrode. The ma¬ 
trix Sin involves phase difference between electron and 
hole converting into two Majorana fermions. This phase 
is denoted as a and is included in the scattering matrix 
as follows 


Vout 


'^in 


^in 


Xout 


L^^nJ 
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'1/V2 

I/V2' 

-gia Q 

_i/V2 


0 e”*" 


(9) 


The fields entering Eq. ( 8 ) are Heisenberg operators at 
a given energy. The Sout matrix was found in Ref. 15 
by means of time reversal transformation of Sin which 
changes the magnetization sign M ^ —M 


Soutm=SU-M). (10) 


We assume here that both Y-junctions have identical ge¬ 
ometries (and, in particular, a is the same). In this work 
we set the value of a as a free parameter which is to be 
found from a particular geometry of the Y-j unction. The 
scattering in the upper Y-j unction in Figure 3 reads as 




Pin 


^out 


7 / 1 + 


Xin 


We match Majorana operators r]in and r]out for the given 
energy 5 as run^e = where the dynamic phase 

= ed/v is accumulated by the Majorana excitation 
in course of the propagation from the lower to the up¬ 
per Y-j unction separated by the distance d. The full 
5'(^S'C=0)-matrix of the left contact is found straight¬ 
forwardly after the exclusion of 77 fields from Eqs. (9), 
(10) and their Hermitian conjugates. To reinstate a non¬ 
zero SC phase of the electrode ^sc (colored as green in 
Figure 3), we employ the transformation -0 
yielding 


charged fermions as linear combinations of uncorrelated 
field operators xi of incident Majorana modes. 

The latter are characterized by the Fermi distribution 
function: 

where 7 ;“^ is the density of states in the yMM channels. 
We assume /cb = 1 everywhere and recover it in final ex¬ 
pressions. Using operator relations, we calculate chiral 
local densities of states and currents ja^jb- The 

linear spectrum of the system (2) means that chiral cur¬ 
rent ji is proportional to charge density and, hence, the 
Josephson current j is given by 

j = 3a - jb = -ev{{lp+lpa} - (14) 


The positive direction of the current is defined from left 
to right SC, bias voltage is zero and the SC phases on the 
electrodes are equal to ±4>/2, as shown in Fig. 2. While 
solving for Dirac T/i^^^-operators, we take into account the 
dynamic and the Aharonov-Bohm (j)AB phases 


= <(.AB = 7 rT. 


The Thouless energy Et here equals to inverse dwell time 
of the interferometer (or the level spacing in the N-region) 


Ey = 


hv 

2L + 2d' 


(15) 


S = 


lA{ke+2o) J_pi "^2^^ kpi{ke-^sc) 


. 2 (xY^sC 




0 


2 ocY^SC 


lpi{ke+^sc) lpi{ke-2a) 

2^ Y2^ 2^ 


( 12 ) 


This unitary S'-matrix has particle-hole symmetry 

0 0 l' 

S{e)=ZS%-s)Z, Z= 0 1 0 

1 0 0 


and acts on Xin,E^ -e)^* describes the par¬ 

tial Andreev reflection in spinless Dirac channel, which 
is combined with neutral Majorana excitations. The An¬ 
dreev part of this process is accompanied by a Cooper 
pair absorption by the SC electrode. 


B. Current-phase relation of transparent junction 

In the following consideration we assume that chiral 
Dirac channels have equal lengths, are parallel to each 
other and separated by the distance d. Using the above 
approach based on S'-matrix (12), we find 7 / 1 -operators of 


Dynamic phase is accumulated by ^-electron or — 5 - 
hole enclosing the interference loop of total length 2L -\- 
2d. Calculations of expressions for are analogous to 
those presented in Appendix B for tunnel junction. 

Within this S'-matrix formalism, valid for energies in 
the range |5| < A, M, we find the following result for the 
C4>R: 


j{^,(f>AB) = ^ j Je{^,4>AB)nF{£,T)de , (16) 

where the distribution function is determined by that of 
the incoming yMMs defined in (13). The spectral current 
Je can be associated with the local densities of states in 
the arms a, 5 via = ev{pi^^ — pi^^). The densities of 
states pi^\ discussed in Section IV C 2, have dimension 
of inverse velocity, hence, is dimensionless. We obtain 

(I>ab) = 

— sin ps sin 4> 

1 + ^ cos((/)AB+4a)+cos ^ ^ _ (^cos(0aB + 4(a) + COS 4>) COS p^ 

(17) 


This is a non-singular 27r-periodic function which re¬ 
flects the structure of the broadened Andreev levels. 
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As we mentioned in the Introduction, conventional 2tt- 
periodicity with respect to $ is related to the assump¬ 
tion of equilibrium (no transport voltage). This means 
that there are no parity effects resulting in anomalous 47r 
Josephson effect. 

The integration in Eq. (16) is reduced to a sum over 
positive Matsubara fermionic frequencies iT{27Tn + tt). 
We note that the dynamical phase takes the values 
cps iT{27Tn + 7t)IEt at these Matsubara frequencies. 
Evaluating the sum, we obtain the C4>R 


j(4>,0ab) = dTT ^ sm4> 

oo ^ 

^ ^ fcgT(l + 2^ • (1^) 

n=o 2e^ — cos{(f)AB + 4(a) — cos 4> 


This is one of the central results of this paper. It de¬ 
scribes the DC Josephson current at arbitrary tempera¬ 
ture and takes into account contributions from the sub¬ 
gap ID states. 

At zero temperature the summation in is 

replaced by an integration over x = and the result 

reads 


jT=o(4>, 0As) 


In I 


A cos(0ab+4q:)+cos $ 

= - 2 -Et —^- 1 -^ 

h co^{(j)AB + 4(a) + cos 4> 


sin4>. (19) 


The C4>R at zero temperature T is presented in Eigure 
4 for the value of the Aharonov-Bohm phase chosen as 
4>ab = —4(a. At this value the first 4>-derivative of the 
current is divergent as oc In |4>| at 4> ^ 0. This divergence 
illustrates the tendency of the C4>R to have spikes at even 
phases 4> = 27rn. In Section IV C we discuss the limit of 
full Andreev reflection, where the C4>R has a sawtooth 
form, also with spikes at 4> = 27rn. 

Two separated Dirac modes connecting the two su¬ 
perconductors form a SQUID loop. In view of the chi¬ 
rality of the junction the Andreev pair belongs to both 
Dirac channels, since a reflection into the same channel 
is forbidden. Considering the junction as a SQUID loop 
controlled by a magnetic flux applied to the — M light 
gray bar in Eig.2, we observe the fractional 24>o = h/e- 
periodic pattern for the critical current. In Eigure 5 we 
plot the critical current as a function of the flux-induced 
Aharonov-Bohm phase (pAB = The curve is sym¬ 

metric under the assumption that Aa = n. (We recall 
that a is related to the geometry of the Y-junctions.) 
The critical current jc in this plot is normalized by its 
maximal value jc,max- At arbitrary a the curve in Eig¬ 
ure 5 would be horizontally shifted. The positions of the 
peaks of jc are given by (jAB.max = -4a + 27rn. 





FIG. 4. T = 0 C4>R of the transparent junction at (^ab = tt. 
Divergent 4>-derivatives at 4> ^ 27rn are observed. 


jc{4^AB) / jc,max 



Aharonov-Bohm phase in N-region, c^ab + 4(a — tt 

FIG. 5. Fractional /i/e-periodic pattern of critical current 
jc{4>AB) at T = 0. The critical current jc(0An) is normal¬ 
ized to the maximal value jc,max and plotted as a function of 
Aharonov-Bohm phase 4>ab = 7r//4>o, where / is the mag¬ 
netic flux. 


IV. TUNNELING REGIME 
A. Scattering matrix of a tunnel contact 

The realization of the tunneling regime of the Joseph¬ 
son junction studied in this paper is presented in Figure 
6. The left and the right tunnel contacts are implemented 
as constrictions of the magnetic material (see Figure 7). 
Similar to the previous Section III, the Dirac channels 
a and b are colored as blue lines and have equal lengths 
la = h = L and geometries. In contrast to the trans¬ 
parent regime, in the tunneling regime the Dirac mode is 
not terminated but rather forms a closed loop. Hence, we 
introduce the total phase acquired over the loop a sum 
of the Berry phase tt and the Aharonov-Bohm phase, 
(j)cxt = tt + 7r//4>o. The constriction of the +M magnet 
plays the role of an insulating barrier. 

We employ the effective ID Hamiltonians (6,7) and add 
a local tunneling term. The resulting Hamiltonian reads 

H = iv j 'tp+d^'tpdx - Y 

tx ^ . (20) 

This Hamiltonian describes counter propagating Majo- 
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rana and Dirac modes which mix at x = 0. The micro¬ 
scopic structure of the contact is accounted for by the 
small tunneling amplitude t v, which is an additional 
parameter of our theory. 

In Appendix A, using Hamiltonian (20), we derive the 
scattering matrix St of the tunnel contact acting on fields 
(V’in.e, Xin,e, ' The Scattering matrix reads 



1 

l+A 


l+A 


1+A^ 


St = 




A ^i^so 
L 1+A^ 


1-A 

l+A 

l+A ^ 


l+A ^ 


l+A 


( 21 ) 


FIG. 7. Structure of the tunneling contact implemented as 
a constriction of the magnetic material marked as dark gray. 
The area of hybridization between the wave functions of the 
counter propagating neutral and charged chiral channels is 
indicated by a bar. 


The formal solution (21) obtained in Appendix A is valid 
for any value of t and the dimensionless parameter A 
characterizing the scattering is given by 

X{t) = ta+ (22) 

x/2v 

Yet, since the Hamiltonian (20) is physically justified 
only in the weak tunneling limit t ^ v, we obtain 
A ^ 1. The amplitude of the domi¬ 

nant process of normal reflection is given by 1/(1 + A). 
The scattering to the Majorana channel scales as r\j vA, 
whereas the Andreev reflection amplitude has the lowest 
amplitude ^ A. 


B. Current-phase relation in the tunneling regime 



FIG. 6. Scheme of the Josephson junction in the tunnelling 
regime. 


Using the scattering matrix St calculated above (21), 
we obtain the following results for the spectral current 
and the C^R (for details see Appendix B): 




—4A^ sin 4> sin 

((1 + A^) cos ifs + cos (pAB — cos 4>)^ + 4A^ sin^ 


(23) 


j{^,(pAB) = 47r sini»y^ 


h 


=0 (1 + A2) cosh (^ "feBn 2 n+l) j 2A sinh (^ "feBn 2 n+l) j + cos (pAB — cos 4> 


(24) 


Here the Thouless energy is given by Et = hv/{2L). 
From the C4>R in (24) we see that at high temperatures, 
T ^ Et^ only n = 0 term contributes to the sum (24). 
In this limit we observe a sinusoidal C4>R and the critical 


current jc is exponentially suppressed: 


i(^) 


AwekBTX^ 


exp{—TrksT/E t) sini>. 


(25) 


fi(l + A)2 
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X 2 

The factor is proportional to the full transparency 

of the junction being the product of the transparencies 
of the left and the right contacts. This result is similar 
to the C^Rs for a conventional S/N/S junction with the 
N-region being a long quantum wire^^“^^. 

In the low temperature regime, T Et^ the C^R is 
no longer sinusoidal and the critical current jc ^ Et de¬ 
cays as 1/L. Flux dependent oscillations of jc are more 
sharper compared to those in the transparent regime. 
The resonant shape of the non-Fraunhofer h/e-periodic 
modulation is shown in Figure 9 for weak tunneling 
A = 0.1. The shift of the maximum of jc by tt results 
from the presence of the Berry phase. 


C. Formal continuation of the tunneling solution to 
the regime of finite transparency 

The solution for the scattering matrix (21) follows for¬ 
mally from the tunneling Hamiltonian (20) at any value 
of t (see Appendix A). Yet, it is only physically justified 
ioY t/v 1, or A ^ 1. We ask ourselves what hap¬ 
pens if we extend (21), formally, to an arbitrary value of 
t/v in the tunneling Hamiltonian (20) and consider the 
scattering matrix (21) at any value of A = tan^(t/v^i;). 
We observe that at A = 1 the C^R obtained for the 
transparent regime (18) is recovered (up to the geom¬ 
etry dependent phase a). In other words, transparent 
contacts formed by Y-junctions correspond to the inter¬ 
mediate strength A = 1 of tunnel contacts, rather than 
to the regime A ^ oc. 

In this subsection we investigate the C^R and the den¬ 
sity of states in the N-region at arbitrary A. In particular, 
we analyze the regime of full transparency, A ^ 00 , where 
the scattering matrix St reaches the unitary limit and 
corresponds to the full Andreev reflection in the Dirac 
channel. 

Note, that the relation between A and t/v, obtained 
in Appendix A, assumes a certain microscopic structure 
of the contacts and, thus, could also be different. We 
discuss this in more details in Appendix A. 


j{^),eET/h 



FIG. 8. Bold curve: C$R found in (24) for tunnel junction at 
dimensionless tunneling strengths A = 0.36. Dashed curves: 
continuation of the results for C4>R (24) to finite transparen¬ 
cies A = 4 and A = 25. The current j(4>) is measured in 
units of the Thouless energy Et with Aharonov-Bohm phase 
4>ab = tt {(jext = 0) and at low temperature T = O.OIFt. 
The A = 25 curve indicates the tendency to the formation of 
spikes at unusual even phases 27rn. 


j c{4^Ab) / j C,max 



FIG. 9. Dimensionless critical current jc{4>AB)/jc,max of 
the tunnel (A = 0.1) and transparent at (A = 1 and 4a = 
tt) junctions as a function of Aharonov-Bohm phase <j)AB = 
27Tfl{hle). Maximum of the critical current is shifted due 
to the presence of Berry phase. Dashed curves correspond to 
results at finite A = 0.5, 2 found from the tunneling approach. 
All figures are plotted at low temperature T = 0.01 Ft. 


reaches the unitary limit 


1. Critical current and C^Rs 


0 0 


As mentioned above, the formal result for C^R for the 
tunneling regime (24) is identical to that for the trans¬ 
parent regime (18) at A = 1, up to geometry depen¬ 
dent phase a. In Figures 8,9 we illustrate the evolution 
of C^Rs and critical currents jc upon increase of the 
transparency parameter A. We plot the results at small 
A = 0.36, related to tunneling approximation, and their 
continuation to higher A = 4 and 25. 

In the limit A ^ 00 the scattering matrix St (21) 


St{\ oc) — 


0 -1 


0 


0 0 


This matrix describes the process of full Andreev reflec¬ 
tion, where an electron converts into a hole with a phase 
shift of wave function equal to the phase of the SC elec¬ 
trode. At T = 0 and A ^ 00 , where jc is independent on 
4>ab^ the C^R (24) shows a sawtooth shape with spikes 








at even phases 

e $ 

jf($) = 4 -^t arctancot — 

= 2^Et{'k - $) for 0 < $ < 27r. (26) 

In conventional ballistic spinful ID S/N/S junctions 
the C^Rs are known to be sawtooth-like^^. In contrast 
to our result (26) the spikes are usually odd it - 1 - 27rn 
and the slope is positive. This difference follows from an 
additional phase of tt acquired by a fermion in a spin¬ 
ful junction after two Andreev reflections from the con¬ 
tacts. For spinful systems, full Andreev reflection is as¬ 
sociated with the time-reversal transformation icjyK of 
the fermion wave function. Being reflected two times 
from the left and right contacts, the spinful fermion gets 
a TT phase. In spinless junction under consideration this 
phase is zero. 


The result for the intermediate transparency A = 1, 
equivalent to the transparent junction, is shown in Fig. 
10 (b). We see that the density of states is strongly 
smeared in this case. 

The continuation of to high transparency (A ^ 1) 
shows that the spectral density is given by 27r-periodic 
narrow lines of Andreev levels. Their structure can be 
found from the singularities of (27) at A ^ oo as 

e;^ = ±£’T(i> + 27rn). (28) 

In Figure 10 (c) we plot the density of states at A = 5 
which consists of half of the full set of smeared Andreev 
levels (28). Note, that level positions are independent 
of the Aharonov-Bohm phase. This follows from the 
fact that the electron and the reflected hole get opposite 
Aharonov-Bohm phases, which compensate each other in 
closed paths. 


2. Density of states 


We calculate density of states at arbitrary A in a- 
branch with the use of ipa from (B 8 ) derived in Appendix 
B. The density of states in chiral Dirac modes are ’halfs’ 
of the spectral current, because = ev{pi^^ ~ if 
coming Major ana fermion modes Xr,i are not correlated. 
We introduce the local retarded Green function (calcu¬ 
lated at an arbitrary coordinate x in the a-branch) as 

(t, t')= - id{t-t'){{lpaix, t), 1 p+ {x, t')}). 


The Fourier transform of this function. 




duo 


uo io 


gives for the density of states 




--IraG, 

IT 




Assuming the density of states of yMM in the SC gap is 
constant, i.e. {{xi-e^Xi.e}) = '^”^5 we obtain 


2A (l + A^ -h cos(0ab - Ps) - cos(^ + ip^)) 

((1 + A^) cos (pe -h COS (^ab — A^ cos <I>)^ + 4A^ sin^ 

(27) 

In Figure 10 we plot the density of states at three 
values of dimensionless tunneling parameter A. At weak 
tunneling, A = 0.3 (see Fig. 10 (a)), the maxima of 
are slightly dependent on 4>. Horizontal lines resemble 
smeared mesoscopic levels of an isolated ID Dirac wire 
of length L. These maxima can be shifted vertically by a 
flux induced Aharonov-Bohm phase because their 

position is given by 

Sn = (27rn -IT - (I)ab)Et. 



FIG. 10. Spectral density of states (27) of the right movers, 
in the upper a-arm of the junction. The energy £ 
is counted in units of Thouless energy Et- Blue and white 
colors correspond to low and high densities respectively, (a) 
weak tunneling regime, A = 0.3; the horizontal lines are rem¬ 
iniscent of the quantized levels of the isolated N-region; (b) 
intermediate transparency, A = 1, the results coincide with 
those for a transparent junction at o = '7r/4; (c) the case 
A = 5 illustrates the continuation to the high transparency 
limit (full Andreev reflection), where bright narrow lines are 
half of full set of smeared Andreev levels. Aharonov-Bohm 
phase <f>AB = —7t/2 for all three plots. 


V. SCATTERING IN TERMS OF MAJORANA 
MODES 

One can describe the scattering between the Dirac and 
Majorana chiral modes by representing the ip operators 
in the N-channel with the help of two auxiliary charged- 
neutral Majorana operators 71 and 72 . These modes do 
not carry charge separately but their superpositions do. 
The different SG phases of the opposite contacts result 
in a fusion of the auxiliary 71 and 72 modes, which is 
responsible for a Gooper pair transfer from one lead to 
the other. 

In this section we discuss the scattering in the Ma¬ 
jorana basis for the tunneling and transparent contacts. 
We analyze in detail the regime A = 1 (keeping a as a 
free parameter), since it corresponds to the case of the 
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transparent junction (18) up to geometry dependent a 
and the Berry phases. 

The Major ana basis is defined through the following 
transformation for both the left and the right contacts 

V’i = (29) 

where the index i here stands for in and out channels. 
The St-msitiix (21) under this transformation at arbi¬ 
trary A reads 




"1 

0 

0 ■ 



7l,OUt 



7l,m 



0 

1-A 

2\/A 



Xout 


l+A 

l+A 


Xin 

Sf2,out. 


_0 

2\/A 

l+A 

1-A 
l+A - 


_72,in •_ 


According to (30) 71 mode is always fully decoupled, 
consistent with the scattering theory of Li, Fleury and 
Biittiker^^. In the weak tunneling limit A ^ 1 the modes 
X and 72 scatter into each other with the amplitude 
^ \/A- In the opposite unitary limit A ^ oc all modes 
are decoupled from each other but both x and 72 invert 
signs meaning that the corresponding Dirac fermions ex¬ 
perience the full Andreev reflection. In the intermediate 
case A = 1 (Figure 11 , left) the modes x and 72 fully 
convert into each other: ^ 2 ,out = Xin and Xout = -72,in- 
The presence of the scattering phase a in the scattering 
matrix ( 12 ) of the transparent Dirac-Majorana contacts 
changes the situations considerably. We apply again the 
transformation (29) to -^-operators and obtain the fol¬ 
lowing S'-matrix in the Majorana basis (setting d = 0) 


71,014^ 


■ cos^ a — sin (+ — ' 


X^^in 

Xout 

= 

— sin (+ 0 — cos a 


Xin 

Sl2,out- 


_ 517 ^ cos a - sin^ a_ 


Sf2,in- 


(31) 

This matrix coincides with that of (30) in the limit A = 1 
only if a = 0. For the other values of a all the modes 
are mixing with each other. At arbitrary a there is mix¬ 
ing between all of the modes except of Xin and Xout- 
Mixing between Xin and Xout is possible if we a add sec¬ 
ond contact with a different SC phase. We illustrate the 
scattering in the particular case oi a = ± 7 r /2 in Figure 
11 (right), where the incoming x-iaode converts into 71 
mode and vise versa. Note that, in this case, 72 mode is 
not converted to the others but gets an inversion of the 
sign, 72 ,owt ~ 'T 2 ,in- 


VI. SUMMARY 

To conclude, we analyzed two limits of ID long ballistic 
Josephson junctions where the leads are formed by gap¬ 
less ID chiral Majorana channels.These junctions can be 


\ 

f Xin 

+M 



XlnOUt 


1 

X2,out 

^SC ‘ 

1 

! 

-M 


\ 

T2,m 

-- 


f 

71,m 


\Xout 




FIG. 11. Scattering in normal Dirac and neutral-charge 

channels x written in terms of auxiliary Majorana modes 
71 , 72 . Left: tunnel contact described by iSt-matrix (21) at 
A = 1. Right: transparent contact described by ^S-matrix 
(12) at o = ± 7 r/ 2 . 


realized as hybrid structures involving 3D topological in¬ 
sulator surface in proximity with s-wave superconducting 
electrodes and magnetic materials. The normal region of 
these setups is formed by two chiral Dirac modes spaced 
by a magnetic material. Such a normal part is effectively 
spinless because the spin textures are locked to profiles 
of magnetic domain walls. 

In the first part of the work we have calculated the 
C^R in the regime of high transparency. In this limit 
the left and the right contacts of the junction consist 
each of two Dirac-Majorana converters built by mag¬ 
net/superconductor interfaces. We find that this sys¬ 
tem has a continuous spectral current, a 27r-periodic non- 
sinusoidal C^R at low temperatures and an h/e-periodic 
dependence of the critical current on the magnetic flux. 
The critical current amplitude at zero temperature is 
given by Thouless energy which is proportional to inverse 
dwell time of the normal region. The junction can act as 
a SQUID because the two Dirac channels are spatially 
separated. An Andreev pair in this case appears to be 
spatially non-local, which offers a possibility of inducing 
an internal Aharonov-Bohm phase. This leads to one of 
our central results: the critical current shows fractional 
h/e-periodic pattern. 

In the second part of the work we have studied an¬ 
other realization of the junction where Dirac and Majo¬ 
rana channels are coupled through a tunnel barrier. The 
C^R of such a tunnel junction was found in terms of 
the tunneling Hamiltonian approach and corresponding 
scattering matrix. The resulting critical current patterns 
show sharper resonant peaks as compared to those for 
the transparent junction. We have also studied the for¬ 
mal extension of the tunneling solution to the high trans¬ 
parency regime and compared it with the one obtained 
in the first part of the work. 

Interferometers involving chiral Majorana modes^^“^^ 
could find their applications as measuring devices of 
topological qubits^^. In this paper we have explored 
a dual setup, in which the equilibrium Josephson cur¬ 
rent is carried by interfering chiral Dirac electrons be¬ 
tween chiral Majorana leads. Embedded into schemes 
with vortices and/or magnetic or SC islands supporting 
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zero-energy Majorana pairs, a Josephson based quantum 
interferometer might be of advantage, e.g., for perform¬ 
ing quantum readout of topological qubits. 
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Appendix A: Derivation of the scattering matrix of 
a tunnel contact 


The local form of tunneling Hamiltonians (20) is an ap¬ 
proximation, because, microscopically, such contacts are 
formed by constrictions of finite size as depicted in Fig¬ 
ure 7. If the constriction is wide and the hybridization 
is negligible, the parallel chiral channels along SC/+M 
and +M/—M junctions have opposite (orthogonal) spin 
textures. In the tunneling area, shown as the constric¬ 
tion of the magnetic material in the black bar, the eigen¬ 
functions can be approximated by superpositions of the 
eigenfunctions of the isolated Dirac and Majorana chan¬ 
nels, ^(^). This means that an electron 

coming into the tunneling area starts to oscillate between 
hole- and Majorana-like states with a rate, proportional 
to a hybridization of the channels, estimated as Me~^^y , 
where dy is a transversal constriction length. Finally, the 
incoming electron scatters into a superposition of outgo¬ 
ing electron, hole and neutral excitation in the Majorana 
channel. From this qualitative picture one can conclude 
that the scattering matrix should be periodic with re¬ 
spect to the phase of these oscillations, given as a product 
of hybridization energy and the constriction dwell time 

Here, we derive this periodic behavior of the St-matiix 
from the local tunneling Hamiltonian (20). We employ 
the Heisenberg equations of motion 


■ {vdo, - dt)ip{x) ■ 


■ -(Pix) ■ 

{vdx - dt)xi-x) 

= it6{x)T 

X(x) 

Xvd^c - dt)i^+{x)_ 


Jp+{x)_ 


where the tunnel matrix is given by 


r = 


^ -i ^sc 
0 e * 2 


i ^^SC „ 

e* 2 0 


0 -e* 2 


-e ^ 2 


(A2) 


We write x(—x) in these equations in order to make the 
chirality of the Majorana mode the same as that of the 
charged channel. In this representation we can consider 
incoming and outgoing states as those at x < 0 and at 
X > 0 respectively. 

It follows from the ^-integration of (A1,A2) around 
the point of contact x = 0 that the relation between the 
tunneling matrix T- and the scattering matrix St reads 

St = exp [i{t/v)T]. (A3) 


Calculating the exponent we obtain 


1 


A p—i^so 

1+A 

1+A ^ 

1+A^ 

1+A ® 

1-A 

1+A 

1+A ® 

A ^i^so 


1 

L 1+A^ 

1+A ^ 

1+A 


where the dimensionless tunneling strength A(t) is given 
by 

A(t) = tan^ —(A5) 
y2v 

This tunneling St-maiiix is unitary and obeys the 
particle-hole symmetry St{s) = ZS^{—e)Z like S. The 
eigenvalues of (A4) are given by ^ 

The solution method leading to (A4) and (A5) is not 
universally accepted. Rather, a different ansatz was used 
in various problems on transport in ID systems such as 
tunneling between edge states of QHE^^, impurity scat¬ 
tering in Luttinger liquid at ^ = 1/2^^ or resonant An¬ 
dreev reflection from zero-mode Majorana bound state^^. 
Eollowing these works we should have taken the local ip- 
and y-operators at point x = 0 as 


pj = 


'^(-0) -h '0(+O) 
2 


x(-0) + x(+0) 
2 


(A6) 


Relations (A6) produce then the same solution for the 
scattering matrix as in (A4). However, the parameter A 
is now different and is given by 


m = 


2^2' 


(A7) 


The two solutions coincide in the weak coupling limit 
t/v 1. Eor larger values of t the difference is substan¬ 
tial. Eor instance, the unitary limit A ^ oc is achieved 
with (A6) and (A7) at t ^ oo. In contrast, with (A3) it 
is reached at tn = j2 -1- iin)v. 

We conjecture that the solution method, leading to 
(A3) applies if the constriction is smooth enough, so that 
the validity of the low energy description provided by (1) 
and (20) is not violated in any point of the constriction. 
Then, Eq. (Al) is solved as a regular differential equation. 
On the other hand, the ansatz (A6) is probably applicable 
for sharp enough constrictions. This ambiguity should be 
resolved by solving 2D BdG equations for the constriction 
geometry. 
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Appendix B: Derivation of C$R in the tunneling 
regime 

In this appendix we present technical details of the 
calculation of the Josephson current. We find the C^R 
as the difference between the chiral currents in a and 
b arms by using the tunnel junction scattering matrix 
St, given by (A4). We take the 1-st and the 3-d lines of 
(A4), disregarding the outgoing Majorana field Xout (line 
2), and obtain 




JPout,—£ 


1 

Ae-‘'*’sc' 


Yin,£ 



l+A 

l+A 


iV^ 

e 2 

Xe'^^sc 

l+A 

1 

l+A _ 


'ib^ 

^in,—E 

+ ^^^X^n,e 

A^SC 

—e^ 2 


(Bl) 

The matrix in the r.h.s. of Eq. (Bl) contains ampli¬ 
tudes of normal and Andreev reflection, while the last 
term describes coupling with the equilibrium Majorana 
in-channel. This term is responsible for the spectral cur¬ 
rent (23) being continuous, due to the gapless spectrum 
of the incoming lead Majorana mode. 

To proceed we rewrite Eq. (Bl) as follows 


X = —1//2. Using (B2) and its Hermitian conjugate at 
5 ^ —5 as well as the property of the real Majorana field, 
xt,-e =xi,s, we derive 


7o 
— A/f* 


-M^I2 

7o 


^s,-L/2 


iV^V^/2Xi,e- (B3) 


Here 70 is the 2x2 identity matrix in the left upper 

. <^ . <^ • ^ -I rT-j 

block and 7^/2 = [e'^^ ^ . Note, that 

for the left contact 6-components of in-fields 

while a-components are oni-fields: Ya,-Lj 2 = Yout and 
'06,-L/2 = Yin- The rank of the 4x4 matrix in (B3) 
equals 2 and the eigenvalues are equal to 0, 0, 2, 2 regard¬ 
less of the values of A and 4>. The condition for the right 
contact is obtained from (B3) by the replacements 


^ ^ Xl Xr, ^e-LI2 e,LI2- (B4) 

The Pauli matrix introduced in (B4), mixes the 2x2 
blocks in (B3) acting as a particle-hole transformation, 
while Xx acts inside of the blocks. The product s^Xx 
interchanges a and b indices in '^s,x- We find that the 
fields in the middle of N-region, = '^e,x=o^ can be 
expressed via the operators at the ends of the Dirac chan¬ 
nels, Tg _/^/2 and X./25 using the dynamical and geo¬ 
metrical phases: 


Yout,i 


= 


Yout,—E 


iX^2XXin,e 


The matrix in (B2) is given by 

'Xei-i’sc i_x 


= 


1 + A -Ae-*‘*’sc 


^SC 

e 2 


e 2 




£,-L/2 




e,0: 


^E,L/2 — ^x'yx'^E^^xXx'^ £,0- 


(B5) 

(B6) 


72) Here = diag[e-*^'/^, and I' = 

diag[e“*'^'**/^, The external 

phase Yext equals to the sum of Aharonov-Bohm and 
Berry phases and the dynamical phase accumulated 
an electron of energy 5 or a hole of energy — 5 , enclosing 
the interference loop of the length 2L. Here 


Using (B2) we formulate the boundary conditions for 
the left and right contacts. We introduce the field = 
[Ya,E^Yt-£^ Ya -E^ Yh^E^^i which depends on 5 and on the 
coordinate x along the Dirac channels. Indices a and b 
stand for the upper and lower Dirac modes (see Eigure 
6). We start from the left contact where 4>5'c' = ^/2 and 



hv 


(B7) 


The relations (B3, B4) together with (B5, B6) make the 
problem of finding the four components of ^^(0) alge¬ 
braically close. The result for the first component reads 


xpa,s = V^e" 


: ^ + ‘P£+0ea;t 




^ecct 


sm 


- Asin^^ 


^) Xi,s + (sin ■ 


. j 4>ext+^ . 

Ae* 2 sm XX 


Xr,E 


(1 + A^) cos — cos Yext — COS 4> + 2iX sin 


(B8) 


The Dirac field in the 6-channel is obtained using the 
geometrical symmetry of the setup and is given by 

Xr, Xl)='>Pa,e{-^, Xl, Xr)- (B9) 

The operator relations (B8) and (B9) are used to calcu¬ 


late 

j = 3a - Jb = -ev j {pPt,ei’a,e) “ (^eV’b.e)) ds, 
which results in the C4>R presented in (24). 
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